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Abstract 



We show that a dense spectrum of chaotic multiply-excited eigenstates can 
play a major role in collision processes involving many-electron multicharged 
ions. A statistical theory based on chaotic properties of the eigenstates enables 
one to obtain relevant energy-averaged cross sections in terms of sums over 
single-electron orbitals. Our calculation of the low-energy electron recombi- 
nation of Au 25+ shows that the resonant process is 200 times more intense 
than direct radiative recombination, which explains the recent experimental 
results of Hoffknecht et al. [J. Phys. B 31, 2415 (1998)]. 
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In this paper we give a quantitative explanation of the puzzle of electron recombination 
with Au 25+ . We also demonstrate how to calculate the contribution of "chaotic" multiply- 
excited states of the compound ion, which mediate electron recombination with complex 
many-electron ions. 

Experimentally, this process was studied recently at the UNILAC heavy ion accelerator 
facility of the GSI in Darmstadt |l| . In spite of a high energy resolution p| the measured 
recombination rate did not reveal any resonances and only showed two broad structures 
around 30 and 80 eV. However, its magnitude at low electron energies e ~ 1 eV exceeded 
the radiative recombination (RR) rate by a factor of 150 although the observed energy 
dependence was close to to that of RR. 

It is well-known that recombination rate can be enhanced by dielectronic recombina- 
tion (DR). In this process the incident electron is captured in a doubly-excited state of the 
compound ion, which is then stabilised by photoemission. Suggested originally by J. Say- 
ers, it was first considered by Massey and Bates [|| in the problem of ionospheric oxygen. 
Later DR was found to be important for the ionization balance in the solar corona and 
high-temperature plasmas on the whole ||. Electron- ion recombination has been measured 
directly in the laboratory since early 1980's ||. More recently the use of heavy-ion accel- 
erators and electron coolers of ion storage rings has greatly advanced the experiment J7|,Q . 
Recombination rates for many ions have been measured from threshold to hundreds of eV 
electron energies with a fraction-of-eV resolution |]9|-[15| . For few-electron ions the measured 
rates were in good agreement with theory which added the contribution of DR resonances 
to the RR background, e.g., He + Hl6| , Li-like C 4H 
T7| . However, more complicated ions, e.g., Au 514 



15] and Ar 15H 
111 and U 28+ 



MM 



and B-like Ar 13+ 



12 



showed complicated 

resonance spectra and strongly-enhanced recombination rates at low electron energies. In 
particular, in U 28+ the theory was able to explain the main resonant features in the range 
80-180 eV, but failed to identify the resonances and reproduce the rate at smaller energies 
|18| |. The situation in Au 25+ |TJ looks even more puzzling. 

In Ref. [19] we suggested that electron recombination with Au 25+ is mediated by complex 
multiply- excited states of Au 24+ , rather than "simple" dielectronic resonances. Electrons 
could be captured in these states due to a strong configuration interaction in this open-shell 
system (the ground state of Au 24+ is 4/ 9 ). The single-particle spectrum of Au 24+ does not 
have large gaps, see Fig. [I]. Using the single-particle orbitals we generated many-electron 
configurations, evaluated their energies and estimated the energy density of multiply-excited 
states ||19||. Owing to the "gapless" single-particle spectrum, the density increases rapidly 



as a function of energy, as described by the Fermi-gas-model ansatz 



p{E) = AE- u exp(aVE) 



(1) 



with A = 31.6, 



v 



1.56, 



3.35 a.u. 21 



where E is the energy above the ground state 
in atomic units used throughout the paper. 

The excited states are characterised by their total angular momentum and parity J 77 , 
and are 2 J + 1 times degenerate. Therefore the total level density is a sum of partial level 
densities: p(E) = J2j*(2J + ^-)pj 7t (E). The spectrum of Au 24+ near the ionization threshold, 
E = I ~ 27.5 a.u., contains J from ^ to ^ EH- Their distribution is in agreement with 
statistical theory, p^ oc /(J), where 
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f(J) 



2(2J+ 1) 
(2J m + 1) : 



exp 



(2J+1) 2 
"2(2J m + l) 



(2) 



J m being the most abundant value of J p0| , p2| . Numerically we find J m ~ |. Using Eq. (g) 
we can estimate the partial densities by pj* = f(J)p/(2J + 1), where (2J+ 1) is an average 
over f(J). For the most abundant J values this leads to pj^{E) = Aj^E~ v exp(a\/E~) with 
Ajtv « 0.15. Near the ionization threshold we have pj* ~ 3.6 x 10 4 au, which means that the 
spacing between the multiply-excited states with a given J n is very small: D = 1/ pj^ ~ 1 
meV. This would explain why autoionising resonances could not be observed in electron 
recombination with Au 25+ ]nj. However, the large density is only a "kinematic" reason 
behind the experimental finding, because we have not proved that the electron can actually 
be captured in the multiply excited states. In what follows we analyse the dynamics of 
electron capture and show that the Coulomb interaction between the electrons makes the 
capture efficient and accounts for the observed enhanced recombination rate. 

Taking into account this interaction is the key problem in many-electron processes. In 
general, this can be achieved by constructing a basis of many-electron states from single- 
particle orbitals, and solving the eigenproblem for the Hamiltonian matrix H ik = (<&i\H\<&k), 
which yields the eigenvalues E v and eigenstates \^ v ) = J2k^k\^k) of the system (config- 
uration interaction method). For open-shell systems with a few valence electrons, e.g., 
rare-earth atoms, this becomes an increasingly difficult task. The density of states grows 
very rapidly with the excitation energy, and finding the eigenstates requires diagonalisation 
of ever greater matrices. 

On the other hand, when the level density is high and the two-body interaction is suf- 
ficiently strong the system is driven into a regime of many-body quantum chaos, where the 
effect of configuration mixing can be described statistically. In this case each eigenstate con- 
tains a large number N of principal components C k ~ 1/y/N, corresponding to the basis 
states <3>fc which are strongly mixed together. This strong mixing takes place in a certain 
energy range \E k — E u \ < r spr , where E k = H kk is the mean energy of the basis state and 
r spr is the so-called spreading width. More precisely, the mean-squared value of as a 
function of E k — E u , is described by a Breit-Wigner formula 
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N~ 
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(E k — E v 



) 2 + r s 2 pr /4' 
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with iV = ttT spt /2D fixed by normalisation J2h 



a 



/ 



c 



dE k /D = 1. The decrease 



E k — E v \ > r spr is a manifestation of perturbation theory: the admixture of 



of Cf } for 

distant basis states is suppressed by large energy denominators. Apart from this systematic 
variation the components behave as Gaussian random variables. 

This picture of many-body quantum chaos is supported by numerical studies of nuclei 
|24| , complex atoms |25] and ions jn| . In particular, Ref . |19j provides an estimate of the 
spreading width in Au 24+ : r spr w 0.5 a.u. Hence, a typical eigenstate near the ionization 
threshold contains JV~2 x 10 4 principal components. 

For low-energy electrons the contribution of multiply-excited autoionising states (reso- 



nances) to the recombination cross section is 26 
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7r 2 7 + 1 rMr( a ) 



P V2(2J, + 1) (£-£,) 2 + r2/4 5 

where e = k 2 /2 is the electron energy, Jj is the angular momentum of the initial (ground) 
target state, J are the angular momenta of the resonances, e v = E v — I is the position of 
the iAh resonance relative to the ionization threshold of the compound (final-state) ion, and 
rw, rw, and r„ = r<r>+r< a > are its autoionisation, radiative, and total widths, respectively 



| 27| . When the resonance spectrum is dense, cr r can be averaged over an energy interval Ae, 



D « Ae « £, yielding 

2tt 2 2j+i / rwrw \ 

* r k 2 ^2{2J 1 + 1)d\ V ^ + T ^I' [ ) 

where (. . .) means averaging. The fluorescence yield Uf = /(Ty+T^>), fluctuates weakly 
from resonance to resonance (see below), which allows one to write a r = cf c (jf, where 



7v'' ^ (2J+ i)rw 



is the energy-averaged capture cross section, and is the average autoionisation width. 

Unlike complex multiply-excited states ty u , the initial state of the recombination process 
is simple. It describes an electron with the energy e incident on the ground state of the 
target (e.g., Au 25+ 4/ 8 , Jj = 6). The autoionisation width is then given by perturbation 
theory as 

= 2ir\(y u \V\$ i ;e)\ 2 (7) 
= 27vY.ct > C^\^e\V\^ kl )^ k \V\^ i -e) ) 

k,k' 

where V is the electron Coulomb interaction, and the continuum states e are normalised to 
unit energy interval. Averaging over the chaotic states v with E u »s I + e, we make use 
of the fact that their components are random and uncorrelated, which leads to 



r^ = 27r]r|cf )2 



'\($ k \V\$ i; e)\ 2 . (8) 



Being a two-body operator, V can move only two electrons at a time. A nonzero contribution 
to r( a ) is given by the basis states which differ from the initial state e) by the positions 
of two electrons. Therefore, in Eq. (§) we only need to sum only over doubly- excited basis 
states With the help of Eq. (|J) the capture cross section (][]) becomes 

- 7T 2J + 1 r spr 27r|($ d |V^;e)| 2 

ac k 2 % 2(2 J, + l)(E d -I- e) 2 + r 2 pr /4' [ } 

This form makes it clear [cf. Eq. (|4])] that the two-electron excitations $d play the role of 
doorway states for the electron capture process. Since these states are not the eigenstates 
of the system they have a finite energy width r spr . 
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The wave function of a doorway state can be constructed using the creation-annihilation 
operators, = a^a^a 7 |$j), where a = n a l a j a m a and (3 = nplpjpmp are excited single- 
electron states, and 7 = n 7 / 7 j 7 m 7 corresponds to a hole in the target ground state. (Note 
that we are using relativistic Dirac-Fock orbitals nljm.) Of course, to form doorway states 
with a given J the angular momenta of the electrons and ionic residue must be coupled 
into the total angular momentum J. However, the 2 J + 1 factor and summation over J 
implied in Eq. (^) account for all possible couplings, and we can simply sum over the 
single-electron excited states a, [3 and hole states 7, as well as the partial waves Ij of the 
continuous-spectrum electron e. As a result, we have 



Or 



TV ^ (a,p\\V x \\~f,elj) 



X 



(«,/3||^||7,^')-(2A+l)E(-ir +A+1 v \(a,P\\Vy\\elj,j} 



A' 




•V jfs h . 



(10) 



where e a , Bp and e 7 are the orbital energies, and (a, /3||Va||7, slj) is the reduced Coulomb 



matrix element f28||. The above equation is directly applicable to targets with closed-shell 



ground states. If the target ground state contains partially occupied orbitals, a factor 

Ua V i-^rl, (11) 




2 h + l\ 2j a + lJ\ 2^ + 1^ 



where n a , rip, and ra 7 are the orbital occupation numbers in the ground state $j, must be 
introduced on the right-hand side of Eq. (|H]). Steps similar to those that lead to Eq. (|10"D 
were used to obtain mean-squared matrix elements of operators between chaotic many-body 

states mm . 



Chaotic nature of the multiply-excited states ^ u can also be employed to estimate their 
radiative widths r^. Electron-photon interaction is described by a single-particle dipole 
operator d. Any excited electron in ^> v may emit a photon, thus leading to radiative stabil- 
isation of this state. The total photoemission rate can be estimated as a weighted sum 
of single-particle rates, 

r^-E#iH*)i 2 (^fi-^^ , (12) 



a,/3 



3c 3 n 11 \ 2^9 + 1 



where iOp a = Bp — e a > 0, (a||<i||/3) is the reduced dipole operator between the orbitals a and 
P, and (. . .) u is the mean occupation number factor. Since \I/„ have large numbers of principal 
components N, their radiative widths display small 1 / y/N fluctuations. This can also be 
seen if one recalls that a chaotic multiply-excited state is coupled by photoemission to many 
lower-lying states, and the total radiative width is the sum of a large number of (strongly 
fluctuating) partial widths. A similar effect is known in compound nucleus resonances in 
low-energy neutron scattering p0[| . 

It is important to compare the radiative and autoionisation widths of chaotic multiply- 
excited states. Equation (|12|) shows that is comparable to the single-particle radiative 
widths. On the other hand, the autoionisation width r^, Eq. (§), is suppressed by a factor 

iV -1 relative to that of a typical dielectronic resonance. A comparison of Eqs. 



0) 2 
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(I) and ([D also shows that is suppressed as D/T spi . Therefore, in systems with dense 
spectra of chaotic multiply-excited states the autoionisation widths are small. Physically 
this happens because the coupling strength of a two-electron doorway state is shared between 
many complex eigenstates. As a result, the radiative width may dominate in the total width 
of the resonances, ^> T^ a \ making their fluorescence yield close to unity. Our numerical 
results for the recombination of Au 25+ presented below, confirm this scenario. 

The resonant recombination cross section should be compared to the direct radiative 
recombination cross section 

32vr Zf , / Zi\ 



3\/3c 3 k 2 \n k 

obtained from the Kramers formula by summing over the principal quantum number of the 
final state [I5|. Here Zj is the ionic charge [Zi = 25 for Au 25+ ), and uq is the principal 



quantum number of the lowest unoccupied ionic orbital. Note that the direct and resonant 
recombination cross sections of Eqs. ( jl3D and (||) have similar energy dependences. There- 
fore, for the purpose of comparing with experiment we can evaluate the cross sections at one 
low energy, say e = 0.5 eV. This energy is much greater than the thermal energy spread of 
the electron beam and the cross sections can be compared directly to the experiment [|TJ. 

Using Eqs. (pi]) and ( [TT]) and summing over the orbitals in Fig. |TJ and electron partial 
waves up to /in/2, we obtain the capture cross section a c = 23 x 1CT 16 cm 2 at e — 0.5 eV. 
A comparison with Eq. @ shows that the sum which contains the autoionisation width, is 

(27 + i)rw 

ji (2Ji+l)D k ' 

Combining this with D = 3 x 10~ 5 a.u. and taking into account that about ten different 
J and two parities contribute to the sum, we obtain ~ 5 x 10~ 7 a.u. On the other 
hand, a numerical calculation of Eq. ( |T2"D gives rW = 3 x 10" 5 a.u. Therefore, r(' r ) 3> r'*' 
and ujf 1. Hence, the resonant recombination cross section is basically equal to the 
capture cross section: o r ~ 23 x 10~ 16 cm 2 . This value is in good agreement with the 
experimental a^ exp - ) = 27 x 10~ 16 cm 2 |l|], and exceeds the direct recombination cross section 
([Hp, a d = 0.12 x 10~ 16 cm 2 , by a factor of 200. 

In summary, we have shown that a dense spectrum of chaotic multiply-excited states can 
play a major role in the dynamics of electron recombination with many-electron multicharged 
ions, and possibly other processes, e.g. charge transfer in collisions of multiply-charged ions 
with neutral atoms. Based on the chaotic nature of these states, we have developed a 
statistical theory which enables one to calculate energy-averaged cross sections for processes 
which go via such resonances. Applied to the recombination of Au 25+ , the theory shows 
that the contribution of resonances exceeds that of direct radiative recombination 200 times, 
which explains the recent experimental findings 

We thank Prof. A. Miiller for providing experimental data in numerical form. 
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FIGURES 




1 2 3 

Orbital rnomentum 



FIG. 1. Energies of the occupied and vacant single-particle orbitals of Au 24+ obtained in the 
Dirac-Fock calculation. 
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